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Model analysis of the very high energy detections of the starburst galaxies M82 
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In the light of recent results detecting the nearest starburst galaxies M 82 and NGC 253 at 
very high energy (VHE), we present a multi- wavelength model succesfully explaining their gamma- 
ray diffuse emission. The detection of M 82 was presented by VERITAS collaboration in the recent 
ICRC, while the integral flux of the fainter NGC253 has just been published by HESS Collaboration. 
Also, in the 2009 Fermi Symposium itself, Fermi data coming from the direction of both galaxies 
was released, thus confirming their detection in the high energy (HE) regime. Slight and reasonable 
variations in the space parameter of already published models can fully account for the HE and VHE 
emission coming from M 82 and NGC 253, while agreeing with previous data detected from radio to 
infrared. We explore these changes and the implications they have for the cosmic-ray distribution 
in these galaxies. 



I. INTRODUCTION 

Starburst galaxies were anticipated as 7— ray 
sources @, 0, provided sufficient instrumental sensi- 
tivity, due to their enhance star formation and super- 
nova (SN) explosion rate, in dense (gas and dust en- 
riched) environments. SN remnants and shock winds 
from massive stars are supposed to accelerate cos- 
mic rays (CR). Due to their collisions with ambi- 
ent nuclei and subsequent ir° decay, very energetic 
7— rays are produced, and those can in turn be de- 
tected both with space-born and ground-based imag- 
ing atmospheric Cherenkov telescopes. 

Recently, the detection of M82 was presented by the 
VERITAS collaboration at the ICRC 2009, while the 
integral flux of the fainter NGC 253 has just been pub- 
lished by the H.E.S.S. collaboration. During the 2009 
Fermi Symposium, the data collected by the Fermi 
space telescope was released and found coming in the 
direction of both galaxies. 

Slight and reasonable variations in the parameter 
space of already published models [2j, |3j can fully ac- 
count for the HE and VHE emission coming from both 
galaxies, while agreeing with previous data detected 
from radio to infrared (IR). We explore these changes 
and some implications they have for the CR distribu- 
tion in these galaxies. 



II. DESCRIPTION OF THE MODEL 

The aim of this and previous studies is to perform a 
multi- wavelength model for the emission coming from 
the central part of the starburst galaxies M82 and 
NGC 253. To model the highest energy emission we 
use the Q—Dif fuse code, already presented in [l| and 
improved in (2| and [|[ . The computation consists on 
solving the diffusion-loss equation for a steady state 



population of both electrons and protons, taking into 
account losses in each particle population, and also 
secondary production. 

One of the major achievements of this study is pre- 
senting an accurate fit to multi-wavelength data. In 
Fig. [l] a multi-frequency spectrum is overplotted to 
previous radio and IR data. Note that the predic- 
tion from the model is the summed contribution of 
synchrotron, free-free and dust emission, and not an 
actual fit to radio or IR data. 



III. GALAXY M 82 

The near, almost edge-on starburst galaxy M82 has 
a low-mass gas content, mostly concentrated in the 
inner 2 kpc, and presents a high luminosity both in 
the far IR and X-ray domain. As part of the M81 
group, M82 shows hints of an encounter with some of 
its members lGyr ago. As a result of tidal forces, it 
harbors a central (300 pc) starburst. 

VHE 7— ray emission coming from M82 was claimed 
by VERITAS Q during this summer, while the HE 
regime was covered by Fermi space telescope from the 
same galaxy. In Fig. fJJ this data is shown together 
with the spectral energy distribution (SED) of the al- 
ready published model b 3] for a range of parameters 
and some specific outputs to better predict observed 
results (see Table H]). A separate contribution is plot- 
ted coming from each 7— ray generator: neutral pion 
(ir°) decay, bremsstrahlung and inverse Compton ra- 
diation. The latter was computed having the cosmic 
microwave background (CMB), far and near IR pho- 
ton densities as targets altogether (see Fig. [IJ. As 
can be seen, ir° decay contribution dominates at VHE 
energies. 
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FIG. 1: Multifrequency spectrum of NGC 253 from radio to infrared. The observational data points correspond to: Carilli 
(triangles, 1996), Elias et al. (circles, 1978), Rieke et al. (asterisks, 1973), Ott et al. (diamonds, 2005) and Telesco et al. 
(squares, 1980) (and references therein). The results from modelling correspond to: synchrotron plus free- free emission 
(dashed), dust emission (dotted) splitted in a cool (blue, T co id = and a warm (purple, T war m 200K) component, 

and the total emission from radio and IR emission (solid). 



IV. GALAXY NGC 253 



V. CONCLUDING REMARKS 



The also near, barred-spiral starburst galaxy NGC 
253 has been deeply studied through the years. Its 
continuum spectrum peaks in the far IR with a high 
luminosity. Its inner (100 pc) region is characterized, 
as well as M82, by starburst activity. 

The integral flux published by H.E.S.S. constraints 
previously predictions for NGC 253 at the VHE 
regime. Moreover, the first results from the Fermi 
telescope add more information to the spectrum of 
the galaxy at lower energies. A set of curves for the 
SED are specifically plotted in Fig. [3] to achieve this 
low flux, exploring uncertainties in the distance to 
this galaxy and subsequent ranges in the magnetic 
field and diffusive timescale. Apart from diffusing 
away during 10 6-7 yrs (see Table HI]) , particles can 
escape the inner starburst convectively, carried away 
by winds (~ 300 km s" 1 ), and through pion colli- 
sions with ambient gas, in even shorter timescales of 
around a few 10 5 yrs. Since the diffusion timescale 
of the particle depends on the energy, the shorter it 
is, the steeper the 7— ray spectrum gets (the higher 
the losses are). The diffusion coefficient associated is 
~ 10 26 - 27 cm" 2 s" 1 at 1 - 10 GeV, compared to the 
~ 10 28 cm" 2 s" 1 value in our Galaxy. 



Our multi-wavelength model explains reasonable 
well both the HE and VHE emission coming from 
the two closest starburst galaxies M82 and NGC 253, 
within a range of explored parameters. Every com- 
ponent of the emission can be tracked to one and the 
same original CR population, and, in turn, result as a 
consequence of all electromagnetic and hadronic chan- 
nels from the primary and subsequently-produced sec- 
ondary particles. 

CR enhancement present in these starburst galax- 
ies is reflected in the high energy density values that 
can be obtained from the steady proton population. 
Above a proton energy of ~ 1500 GeV (corresponding 
to E 7 ~ 250 GeV), the energy density is around 10 
eV cm" 3 for M82 and similar value for NGC 253. 

Now that the VHE regime has been finally achieved 
by ground-based telescopes, and Fermi has proven 
that detections are possible from lower energies ( 100 
MeV), a full picture of 7— ray emission is begining to 
appear more and more clearly. 
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FIG. 2: Energy distribution of the differential gamma-ray fluxes of M82, exploring a range of uncertainties in supernova 
explosion rate and efficiency to inject energy from SN to CR. The shaded green area corresponds to the original model of 
de Cea del Pozo et al. (2009) . Data points and upper limit correspond to both VERITAS (stars) and Fermi (diamonds) 
detections. 
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FIG. 3: Energy distribution of the differential gamma-ray fluxes of NGC 253, exploring the uncertainty in distance, a 
range of timescale diffusion (to) and possible cutoffs in the proton injection spectrum. The original model from Domingo 
& Torres (2005) is also shown for comparison, as well as data points from Fermi detection (diamonds) and the integral 
flux from the H.E.S.S. detection (square) transformed in differential flux (assuming a range of injection spectra). 
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TABLE I: Physical parameters used in the multi-wavelength model of M82, both in 0] and the present study, in which 
values are specified in order to match with the VERITAS detection. The (small) variations explored here are, in any 
case, inside the former predictions of the original model. Numbering divide the list of parameters in: 1) observational 
values, 2) derived from 1), 3) obtained from modelling, and 4) assumed. SB stands for starburst. 







Ho Ppn Hoi Pn7n of ( 9,nflQ i 


VERITAS-Hriupn moHpl 




Distance 


q Q _|_n q in q A/Tr^r- 
<j.c ^random -LW.o systematic iV1 P^ 






Inclination 


77-LQ ° 
/ / IO 






Haaius od 


oUU pc 






Radius Disk 


7 kpc 






rieignt . j 1 3 


90,0, nn 

zuu pc 






ijras iviass od 


Z X 1U M0 t^il 2 J 






L*as Mass Disk 


i \/ i o& A/f /"ui^ 
/ X 1U M0 t,JllJ, 




w 




l.o X 1U M0 \-tl2) 






TT-? T liininACihr 
LLX, J_j LlllllllUol L V 


a v i n 10 r ^ 

■4: A J_1J IsQ 






SN explosion rate 


u.o yr t^u.i yr ) 


u.z yr — u.o yr 




^\I\T ovnl ncmn onormr 
OIN caJJIUdIUII cllGl cv 


ID 51 prff 






SllSJ pti prffir tran <3TPT"~rprl 1~n C;T? 

kJl'l CHLl si y til dllolCl i. CU LU V • 1 V, 


10% 


1(1% Fi% 




Convective vGlocity 


OUU Kill S 






l)i i tit l~p~m nprati i fp 


45 K 






Inn 1 7pr1 tpTn nprati 1 t*p 


1UUUU IV 






T IniTAriTi Hpncitir 

U I111UI 111 LlLIlSl L V O-U 


lOULi/t 






Dust emissivity index 


1.5 






Emission measure 


5 x 10 5 pc cm" 6 




(3) 


Magnetic field 


120 fiG (270 fiG) 


170 fiG — 210 




Proton to electron primary ratio 


50 (30) 






Slope of primary injection spectrum 


2.1 






Maximum energy for primaries 


10 6 GeV 




(4) 


Diffusion coeficient slope 


0.5 






Diffusive timescale 


1-10 Myr 





Acknowledgments 

The work of E. de Cea del Pozo has been made 
We acknowledge support by grants AYA2006-00530 under the auspice of a FPI Fellowship, grant BES- 
and SGR2009-811. 2007-15131. 



[1] Torres, D.F. (2004), ApJ 617, 966. 
[2] Domingo-Santamaria, E. & Torres, D.F. (2005) A&A 
444, 403. 

[3] de Cea del Pozo, E., Torres, D. F. and Rodriguez Mar- 

rero, A. Y. (2009) ApJ 698, 1054. 
[4] Acero et al. (H.E.S.S. coll.) Science Express, 



10.1126/science.ll78826, (2009) 
[5] Benbow et al. (VERITAS coll.), 0656, ICRC Lodz 
(2009). 

[6] Paglione et al. (2009) ApJ 460, 295 
[7] Persic et al. (2008) A&A 486, 143 



eConf C091122 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



5 



TABLE II: Physical parameters used in the multiwavelength model of NGC 253, as presented both in the previous @] 
paper and in the present study, which explores some variations allowed within the model in order to explain H.E.S.S. 
detection. Numbering divide the list of parameters in: 1) observational values, 2) derived from 1), 3) obtained from 
modelling, and 4) assumed. SB stands for starburst. 
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